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Abstract 

We apply geometric techniques to obtain the necessary and sufficient 
conditions on the existence and nonlinear stability of self-gravitating Rie- 
mann ellipsoids having at least two equal axes. 



1 Introduction 

A Riemann ellipsoid is a relative equilibrium for a dynamical model of a rotating 
self-gravitating fluid mass that remains an ellipsoid at all times, and for which 
the fluid velocity is a linear function of the coordinates. This model, first studied 
and formulated by Dirichlet, is nowadays known as Dirichlet's model. The 
linearity assumption on the allowed motions makes the study of these deformable 
bodies very attractive since it implies that their dynamics is governed by a 
system of ordinary differential equations with a finite number of degrees of 
freedom. These bodies are also known as affine-rigid bodies or pseudo-rigid 
bodies. Dirichlet's model can be viewed as a first order model for the study of 
the shape of the Earth, and the study of the stability of its solutions can be 
used in planetary stability research. 

The study of self-gravitating fluid masses has a long history which can be 
traced back to Newton's times with many contributions from Dirichlet, MacLau- 
rin, Jacobi, Dedekind, Riemann, Liapunov, Poincare and Cartan, just to name 
a few. We can distinguish two classical approaches to the study of the stability 
of Riemann ellipsoids. One initiated independently in the latter part of the 
nineteenth century by Poincare [18j and Liapunov |10) who used expansions 
in spherical and ellipsoidal harmonics to study the stability of MacLaurin el- 
lipsoids and Jacobi ellipsoids, not only under Dirichlet's assumptions but also 
under more general conditions (allowing perturbations not preserving the el- 
lipsoidal shape). The other approach occurred in the middle of the twentieth 
century with the works of Chandrasekhar and collaborators who developed the 
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so-called virial method by applying it to the study of the linear stability of Ric- 
mann ellipsoids. These works are collected in the book of Chandrasekhar |3J 
which constitutes a comprehensive survey on the subject and related problems 
with many historical facts on this model. 

In recent times the subject has had the attention of several researchers, in 
particular in what respects to the application of new formulations and meth- 
ods ([24l [25]) to study rotating deformable bodies, not only subject to the 
self-gravitating potential but also for other potentials modelling nuclei (see for 
instance [13]), or elastic bodies (see [5] and [S]). 

Our aim is to use geometric methods not only to obtain a complete cha- 
racterization of the conditions for the existence of Riemann ellipsoids having 
configurations with at least two equal axes but also to obtain the complete 
description of their nonlinear stability. These geometric methods exploit the 
geometry and symmetries of the problem and its Hamiltonian structure. Some 
works using the same philosophical approach to Dirichlet's model are available 
in the literature for studying several aspects of the problem such as in [2TJ| 
which obtains Riemann's theorem as a consequence of the symmetry alone, or 
[8] where some results on the formal stability and bifurcations of MacLaurin 
spheroids are obtained. 

We do not address the problem of the stability for Riemann ellipsoids with 
three distinct axes for which the self-gravitating potential is an elliptic integral. 
In [B] this has been studied employing numerical analysis techniques. 

We view Dirichlet's model as a Hamiltonian system where Hamilton's func- 
tion, h, has the form kinetic plus potential energy and is defined on the cotangent 
bundle, T*SL(3), of the set of all 3x3 matrices of determinant 1. Furthermore, 
h is invariant for the action of G = Z2 x (SO(3) x SO(3)) on the phase space. The 
SO (3) x SO (3) symmetry reflects the existence of two conserved vector quan- 
tities (by Noether's theorem): the angular momentum and circulation. The 
Z2 symmetry reflects the reciprocity theorem of Dedekind: interchanging the 
angular velocity and vorticity vectors one obtains another (physically different) 
solution for Dirichlet's model with the same geometric configuration. 

It is well known that Noether's conserved quantities are organized as the 
components of the momentum map. In the case of relative equilibria with 
configurations having at least two equal axes (symmetric configurations) the 
corresponding momentum map value can be singular. Recently, one of the 
authors developed in |21j a method appropriated for the study of nonlinear 
stability of this kind of relative equilibria. This construction extends the so- 
called reduced energy-momentum method of Simo, Lewis and Marsden [26j to 
the case of singular relative equilibria and so we will refer to the method [21] as 
the singular reduced energy-momentum method. This is the approach we use in 
this work to study the nonlinear stability of Riemann ellipsoids with symmetric 
configurations. 

The preliminary sections of the paper are organized as follow: in Section d] 
we give the geometric formulation of Dirichlet's model, in Section [3] the sin- 
gular reduced energy-momentum method is briefly reviewed and in Section [4] 
we compute the augmented potential energy for symmetric configurations, as a 
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necessary step towards the stability analysis. 

The main results of the paper are in the following two sections. In Section[5l 
Theorem 15. 1[ we give the complete characterization of all the possible Riemann 
ellipsoids with symmetric configurations. We prove that for Dirichlet's model 
the only relative equilibria with configurations having at least two equal axes 
are: the spherical configuration which has zero angular momentum and circula- 
tion, the MacLaurin spheroids which are oblate spheroids rotating around the 
(shortest) symmetry axis and have angular momentum and circulation aligned 
with it, and the transversal spheroids, which have prolate spheroidal configura- 
tions that rotate around an axis, say n, perpendicular to the (longest) symmetry 
axis and have angular momentum and circulation aligned with n. We also prove 
that there are no symmetric relative equilibria for which the angular velocity 
and vorticity are not aligned with a principal axis of the relative equilibrium 
configuration. That is, there are not symmetric configurations which are not of 
type S in Chandrasekhar's terminology. 

In section|5]we apply the singular reduced energy-momentum method to the 
study of the nonlinear stability of the relative equilibria found in the previous 
section. The main results of this section are theorems 16.11 16.21 and 16.31 giving 
respectively necessary and sufficient conditions for the nonlinear stability of the 
spherical equilibrium, MacLaurin spheroids and transversal spheroids. 

Our results on the eccentricity range for the nonlinear G^-stability of the 
MacLaurin spheroids agree with those already obtained by Liapunov and Poincare 
(see Remark 16.11) . In the works of these authors there is no reference to the 
transversal spheroids, however their existence is acknowledged in pag. 143 of 
Chandrasekhar's book [3]. It is not clear for us what are the results obtained by 
Chandrasekhar with respect to the stability of the transversal spheroids, how- 
ever his use of the virial method only gives linear stability. To the best of our 
knowledge, our result on the nonlinear stability of the transversal spheroids is 
new. 

In conclusion, this work presents, from a purely geometric point of view, 
a self-contained and complete study of the nonlinear stability of all symme- 
tric relative equilibria for Dirichlet's model. At the same time, the richness 
of the model helps to clarify the applicability of the singular reduced energy- 
momentum method, providing also a methodology for its application to other 
models like. 

Acknowledgements. We would like to thank to Centro de Analise Matematica, 
Geometria e Sistemas Dinamicos of the 1ST, Lisbon, for the Portuguese trans- 
lation of original Riemann's paper [TH] made by C.E. Harle. The work of ESD 
has been supported by the Fundacao para a Ciencia e a Tecnologia through the 
Program POCI 2010/FEDER. 

2 Geometric formulation of Dirichlet's model 

Let (M, <C •, • 3>) be a Riemannian manifold (the configuration manifold), G a 
Lie group that acts by isometries on M (the symmetry group) and V € C°° (M) a 
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G-invariant function (the potential energy) . With these ingredients we construct 
a symmetric Hamiltonian system on T*M (which is a manifold equipped with 
a natural symplectic structure) as follows: the potential energy V can be lifted 
to T*M with the pullback of the cotangent bundle projection r : T*M — > M. 
We denote this lifted function also by V. The Riemannian metric on M induces 
an inner product on each cotangent fiber T*M , x G M. Then the Hamiltonian 
is defined as 

h(Px) = ~\\Pxf + V(x), Px eT*M. 

The G-action on M induces a cotangent- lifted Hamiltonian action on T*M with 
associated equivariant momentum map J : T* M —f q* defined by 

(J(P«),0 = <P»,M*)> V£e fl , 

where £m is the fundamental vector field on M associated to the generator £, 
defined by 



This momentum map is Ad*-equivariant in the sense that 3(g-p x ) = Ad*-iJ(p x ) 
for every p x G T*M, g G G. 

The Hamiltonian h is G-invariant for this lifted action (this follows from 
the invariance of the metric and of V). Therefore, due to Noether's theorem, 
the components of J are conserved quantities for the Hamiltonian dynamics 
associated to h. The quadruple (M, <C •, • ^g>, G, V) is called a symmetric simple 
mechanical system. 

Let (M, -c v 3>,G, V) be a simple mechanical system with symmetry. A 
relative equilibrium is a point in phase space p x € T*M such that its Hamilto- 
nian orbit lies inside a group orbit for the cotangent-lifted action. This amounts 
to the existence of a generator £ S g such that the evolution of p x is given by 
e'* • . The element £ is called a velocity for the relative equilibrium and is 
defined up to addition of elements in g Px = Lie (G Px ), where G Px is the stabilizer 
of p x under the cotangent-lifted action. A useful criterion for finding relative 
equilibria in simple mechanical systems is given by the following theorem: 

Theorem 2.1 (Marsden [TT]). A point p x 6 T* M of a symmetric simple me- 
chanical system (M, -C •, • S>, G, V) is a relative equilibrium with velocity £ G q 
if and only if the following conditions are verified: 

1. p x =< £ M (x), ■ >• 

2. x is a critical point ofV^ := V — !(•)(£)), 

w/iere I:Mxg-^g* is defined by (£,I(x)(r])) =<C £a/(^), Vm(x) Moreover, 
the momentum fi — J(p x ) G 0* of the relative equilibrium is given by /.t = I(a;)(£). 

Note that, in virtue of the above theorem, any relative equilibrium is char- 
acterized by a configuration- velocity pair (x,£) G M x g satisfying dV^(x) = 0. 
The map I is called the locked inertia tensor, while the function is called the 
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augmented potential. We indicate for later use that the kernel of I is precisely 
q x , the Lie algebra of G x , the stabilizer of x for the G-action on M. The knowl- 
edge of the pair allows us to compute the stabilizer of the corresponding 
relative equilibrium p x —<t^. 1;m{%), ■ 3> with the formula (see |22j): 

G Plc ={geG x : Ad 3 £ - £ g 0*}. 

Dirichlct's model is a symmetric simple mechanical system for the motion 
of a homogenous and incompressible fluid mass of density p having as reference 
configuration the unit ball centered at the origin in R 3 and subject to the self- 
gravitating potential. The only allowed configurations for this model are linear 
embeddings of the reference ball into R 3 preserving volume and orientation. 
The configuration manifold M for a self-gravitating fluid mass under Dirich- 
let's conditions is then SL(3), the group of all 3 x 3 matrices with determinant 
equal to 1, which is equivalent to the space of orientation and volume preserving 
linear automorphisms of R 3 . In what follows we review the geometric formula- 
tion of Dirichlct's model as a symmetric simple mechanical system on GL + (3), 
the group of all 3 x 3 matrices with positive determinant, with a symmetric 
holonomic constraint. 

The singular value decomposition of any linear automorphism of R 3 allows 
to decompose (non- uniquely) any matrix F G GL + (3) as 

F = LAR T 

where L,R G SO(3), and A is a diagonal matrix with positive entries called 
singular values (the square roots of the eigenvalues of C = F T F). It follows 
from this decomposition that the reference unit ball is mapped by F into a 
solid ellipsoid of equation X • G _1 X = 1, X G R 3 , having principal axes half- 
lengths equal to the entries of A. Physically the matrix L describes the rigid 
rotation of the body in space relative to an inertial frame and R is related to 
the rigid internal motion of the fluid with respect to a moving frame. Then 
A is an orientation preserving dilatation of the original reference body into an 
ellipsoid with principal axes aligned with the eigenvectors of A. The condition 
on the volume preservation of the total embedding corresponds to impose the 
holonomic constraint det-F = 1 (or equivalently detA = 1), which in turn 
amounts to consider our system as defined on SL(3). 

The tangent space at F £ GL + (3) is isomorphic to L(3), the vector space of 
3x3 matrices. We can define a Riemannian metric on GL + (3) as: 

< 5F 1 ,6F 2 >= T tr{5F?6F 2 ) (1) 

for SFi, 5F 2 G 7VGL + (3), and T is a constant depending on the density of the 
reference body and other physical parameters of the system. In the case of 
interest here, the reference body is a homogeneous unit ball of constant density 
p, and T in (P) is 
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The symmetry group G of our model is the semi-direct product G = Z 2 X 
(SO(3) x SO(3)), where Z 2 = {e, cr}. Several actions of G of interest in this 
paper are: 

(1) The G-action on G: If (7; g, h), (7'; g', h') e Z 2 x (SO(3) x SO(3)) then 

(7; 9, h) ■ (7'; 9, h') = (77'; (g, h) ■ (7 • (5', /i'))> 

where for the nontrivial element cr € Z2, cr - (g', h') = (h',g'), and SO(3) x 
SO (3) acts on itself by the direct product of left matrix multiplications. 

(2) The adjoint representation of G: We identify the Lie algebra of 5*0(3), 
the set so (3) of skew-symmetric 3x3 matrices, with R 3 via the usual 
isomorphism^: R 3 — * so (3): 







' 


-V3 


V 2 


v = (vi,v 2 ,v 3 ) H 


■*v = 


^3 





-Vl 






-v 2 








This is an isomorphism of Lie algebras, i.e (so(3), [• , •]) is isomorphic by ^ 
to (M 3 , x), where [, ] denotes the commutator of matrices and x denotes 
the vector product of vectors in M 3 . 

The Lie algebra of G is then g = M 3 © R 3 and the adjoint action is given 
by 

Ad (7;g!/l) (£ L ,6?) =nr-(g'€L,h-£ R ) 

where a ■ (£l,£r) = (£rj£l) and g ■ is the rotation of £l by g (and 
similarly for 

Using the standard inner product in R 3 we also identify g* with R 3 © R 3 . 
Under this identification it follows easily that the coadjoint representation 
has the expression 

Ad ( 7 ; ff ,fc)-i(Mi>Mfl) =J-{g-HL,h-n R ) 

(3) The G-action on GL+(3): 

(e; L,R) ■ F = LFR T , (cr; L, R) ■ F = RF T ' L T , 

Note that the Z2 transposition symmetry on SL(3) (first noticed by Dedekind) 
maps a rigidly rotating configuration without internal motion into one that is 
stationary in space but with rigid fluid internal motions. That is, for a given 
ellipsoid there is an adjoint one, obtained by transposition. These adjoint type 
of ellipsoids are called Dedekind ellipsoids. More generally, the transposition 
symmetry interchanges external rotations and internal motions for any solution 
of Dirichlct's model. 

Any G-invariant function / on GL + (3) can be written as 

f(F) = f(I 1 (F),I 2 (F),I 3 (F)), 
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where / : M 3 — * R, and I\ , 1% and I3 are the three principal invariants of a 3 x 3 
matrix, given by 

h(F) = tr(S) (3) 
12(F) = i(tr 2 (S)-tr(S 2 )) (4) 
I 3 (F) = det(S) (5) 

with S = FF T (also valid interchanging S with C = F T F). Note that I\,h,h 
are G-invariant and that 1 is a regular value of I3 . Hence we haveSL(3) = 1^(1) 
as a G-invariant submanifold of GL + (3). The condition I3 = 1 is the holonomic 
constraint of the model. Note also that the restriction of a G-invariant function 
/ £ G G (GL+(3)) to SL(3) is given by 

f(F) = f(h(F),I 2 (F),l), 

for F £ SL(3) C GL + (3). Therefore any G-invariant function on SL(3) can be 
written as h(F) = h{I\(F),h(F)), with h : M 2 — > R. Any G-invariant function 
on SL(3) can be then extended invariantly to GL + (3) by declaring it to be 
independent of I3 (F) . From now on we will drop the tildes and hats from the 
corresponding functions unless there is risk of confusion. 

Since G acts on GL + (3) by isometries with respect to ([1) , the induced metric 
on SL(3) (which we will denote by the same symbol) is also G-invariant. For 
later use we recall that tangent vectors to SL(3) can be seen as tangent vectors 
to GL + (3) satisfying the linearization of the constraint I3 = 1. In other words, 

7>SL(3) = {SF £ L(3) : tr(F -1 «yF) = 0}. 

Dirichlet's model is governed by a Hamiltonian function of the form kinetic 
plus potential energy on the phase space V = T*SL(3) given by 

h(pF) = ^\\ P F\\ 2 + V(F), p F £T£SL(3). 

Here \\pf\\ is the norm of the covector pp (seen as a 3 x 3 matrix) relative 
to the metric on SL(3). The potential energy V for a self-gravitating body of 
homogeneous density p under Dirichlet's assumptions is given by restricting the 
function 

where F £ SL(3), R = j^n 2 Gp 2 , G is the gravitational constant and 

A(F) = Vs 3 + /i(^)s 2 + / 2 (F).s + l. (7) 

The quadruple (SL(3),< v >,Z 2 x (SO(3) x SO(3)),V) defines a symmetric 
simple mechanical system on SL(3). 
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The infinitesimal generator for the G-action on M = GL + (3) (and on SL(3)) 
corresponding to £ = (£lj£r) € M. 3 X M 3 is: 

M^) =i| t = (expCexp^)- J F = &F - F&- (8) 

The vectors £l and £r are respectively the angular velocity and vorticity of the 
fluid motion. We denote the momentum value of pp by J(pf) = (j , c) . The 
components j and c are respectively the angular momentum and circulation of 
the instantaneous state pf (see for instance [2U] or [5]). 

A Riemann ellipsoid (a.k.a. an ellipsoidal figure of equilibrium) is a solution 
of the Hamiltonian system defined by Dirichlet's model with angular velocity, 
vorticity and principal axes lengths all constant. In our setting, Riemman el- 
lipsoids correspond exactly to relative equilibria of the underlying symmetric 
simple mechanical system. Therefore, a Riemann ellipsoid is represented by 
a triple (F, £l , £,r) , where F € SL(3) is the configuration matrix and the Lie 
algebra element (£l,£r) € is the angular velocity- vorticity pair. 

3 The singular reduced energy- momentum method 

In the last years, several are the works studying the stability of relative equi- 
libria of Hamiltonian systems by exploiting their symmetry and the geometric 
properties of their phase space (see for instance [2J, [T5], and [TT] for a overview). 
Most of these methods can be used to test the stability of relative equilibria ly- 
ing in singular level sets of the momentum map, for instance [7] and |14| under 
the hypotesis of G^ compact and [13], [16] using topologic properties. This 
observation is important since, as we will see, the class of symmetric Riemann 
ellipsoids known as MacLaurin spheroids corresponds precisely to non-trivial 
relative equilibria for Dirichlet's model having singular momentum values. We 
refer the reader to [16) for a comparison of the applicability of several existing 
methods. 

The generally adopted notion of stability for relative equilibria of symmetric 
Hamiltonian systems is that of G^-stability, introduced in [15] and that we now 
review in the context of symmetric simple mechanical systems. This notion is 
closely related to the Liapunov stability of the induced Hamiltonian system on 
the reduced phase space. 

Definition 3.1. Let (M, <C •, • ^>,G,V) be a symmetric simple mechanical 
system and p x £T*M a relative equilibrium with momentum value fi = J(p x ). 
We say that p x is G M stable if for every G ^-invariant neighborhood U C T* M of 
the orbit G M • p x there exists a neighborhood O of p x such that the Hamiltonian 
evolution of O lies in U for all time. 

In the early 1990's ([2B]) a tool known as the reduced energy- momentum 
method has been developed, providing sufficient conditions for the stability of 
relative equilibria of a simple mechanical system under the hypothesis that its 
momentum is a regular value of the momentum map. This method is espe- 
cially well-suited for simple mechanical systems since in incorporates all of their 
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distinguishing characteristics with respect to general Hamiltonian systems. Re- 
cently, based on the characterization [17] of the so-called symplectic normal 
space N for a cotangent-lifted action, the reduced energy-momentum method 
was generalized in 21J to cover also the case of singular momentum values. 

In this section we outline the implementation of this singular reduced energy- 
momentum method following [21] . Our setup will be as in Definition 13.11 and 
Section 2. In particular we will fix a relative equilibrium p x with configuration- 
velocity pair (x,£) and momentum /i. We will also assume that the G-action on 
M is proper and that there exists a G M -invariant complement to g^ in g. These 
last two conditions are always satisfied for any relative equilibrium in Dirichlet's 
model due to the compactness of G. 

We start by stating some key observations: First, by equivariance of J : 
T*M -> g* and r : T*M -> M, one has G Px C G x and G Px C G M . In fact, it is 
not difficult to prove the characterization 

G Pm =G x nG ll . (9) 

We remark that the above formula is not valid in general for covectors p x other 
than relative equilibria. 

Second, also by equivariance of r together with the Bifurcation Lemma (see 
[I]), if is a singular momentum value then q x ^= {0}, in which case \i € 
(Qx)°- Third, the properness of the G-action implies that G Vx is compact. This, 
together with ([9]) allows to define the following G Px -invariant splittings: 

0fi = 0p* © P and = Sx © p ffi t, (10) 
which by duality induce the splittings 

*t=Q* P * ®p* and 0* = 0;©p*©f*. (ii) 

Here p and t must be chosen in such a way that l(x)(p) C t°. Note, from the 
definition of the locked inertia tensor in Theorem 12. 11 that kerl(x) — g x , so the 
restriction 

I o =I(.x)| pel :p©t->( 0;E ) o =p*©t* 

is a Gp x -equivariant isomorphism. Then the condition on the above splitting 
is that p and t must be orthogonal with respect to the inner product on p © t 
induced by I . 

We will denote generically the linear projections associated to the splittings 
(p~0|) and (TrTj) by the letter P with an appropriate subindex. For instance P p : 
g — > p or P t » : g* — > t*. It is a consequence of Noether's theorem that £ G g M 
and so we will denote by £ = Pp(£) € p (the orthogonal velocity of the relative 
equilibrium) . 

In this work we only need a particular version of the singular reduced energy- 
momentum method. For, consider the following definitions: 

(1) Let S be the linear orthogonal slice for the G-action on M at x, i.e. 

S = (T X {G ■ x)) x e T X M. 
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(2) Define the subspace q M C g as 

^ = {7 G t : P sJ [ad>] = 0} . 

(3) Define the space of internal variations 

Smt = {7m +a: 7°£^«eS, and (DI • (F) + a)) (£ x ) G p*} . 

(12) 

(4) For any Vi,v% G T X M, the correction term is the bilinear form on T X M 
defined by 



corr 1*) = (Pp-et- [(DI • vi) (£)] , V (Vet- [(DI • «a) (£)]) 



(5) The Arnold form Ar : q^ 1 x q^ 



is defined by: 



Ax (71, 72) = (ad*^, I (ad 72 /i) + P pe t 



ad-. 



(Vm)]) 



(13) 



(14) 



The following theorem (Corollary 6.2 of [5T]) is the synthesis of the singular 
reduced energy-momentum method. 

Theorem 3.1. Let p x G T*M be a relative equilibrium with configuration- 
velocity pair (x, £) G Mxg and momentum fi G 0* suc/i £/iai dim(G--F) < dimM. 
Lei = P p (£) fee ite orthogonal velocity. If the Arnold form is non-degenerate 



at p x and (d x V^± + corr t± (x) J 
riwm is G ^-stable. 



is positive definite, then the relative equilib- 



When the Arnold form is non-degenerate it is also shown in [21] that the 
symplectic matrix of the symplectic normal space at p x has a particularly simple 
block-diagonal expression. We quote this result which will be essential in the 
proof of the linear unstability of some Riemann ellipsoids. 

Theorem 3.2. If the Arnold form is non-degenerate at the relative equilibrium 
p x with configuration-velocity pair (x, £) and momentum fj,, the symplectic nor- 
mal space at p x , is symplectomorphic to N = q M © Si nt © S* equipped with the 
symplectic matrix 



q" 




S* 









* T 




1 





1 






where 



2(71,72) = -(/x,ad 7l 7 2 ), 



t>{lAlM{x)+b)) = {pLM 1 l b ) 



and is the one-form defined by x^ (A) =<C £,m,X for all X G X(M). 
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Remark 3.1. The G-invariance of V and <C •, • ~3> imply the following property. 
If is a relative equilibrium, then the orbit {g ■ x, Ad 9 £) for every g 6 G 

consists of relative equilibria with the same stability or unstability properties. 

Remark 3.2. The reason why in the previous section we looked at Dirichlet's 
model as a simple mechanical system holonomicaly constrained is that the un- 
constrained space GL + (3) is an open domain of the vector space L(3), and then 
the implementation of the reduced energy-momentum method is easier than if 
one is working directly on SL(3). In view of the survey of the method, the 
strategy will be to use the trivial extension of the self-gravitating potential to 
GL + (3) and consider its augmented potential with respect to the locked inertia 
tensor corresponding to the original Riemannian metric on GL + (3). Then we 
further augment this augmented potential with the constraint function /3 and 
Lagrange multiplier A. Denoting by 

the resulting twice augmented potential, we have: 

(1) Relative equilibria for Dirichlet's model correspond to triples (F,£l,£r) 
with F € GL + (3) , (£l , £r) G K 3 x R 3 such that the following two equations 
hold: 

dV^ R (F) = 0, and det(F) = 1. (15) 

(2) The stability test now follows from the following observation. If we call 

^int ^ an d ^fn't*' 3 '' tne spaces of internal variations for GL + (3) and SL(3) 
associated to the triple (F, £l,£r), according to (|12p. we notice that 

^ (3) =Sp n f (3) nkerT F / 3 . 

Therefore, according to the general method, and standard Lagrange mul- 
tiplier theory, to conclude stability it suffices to study the eigenvalues of 
the bilinear form 



(d F V^ R) x+corr ( ^ R) x(F) 



,SL(3)J 



where sf^ 3 ^ is seen as a vector subspace of From now on we 

will omit the superindex SL(3) for the space of internal variations. 

4 The augmented self-gravitating potential for 
symmetric configurations 

In this section we compute the augmented potential 

V iL , 6R =V-^ L ,CR),I(F)(^^R)) 
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in the unconstrained configuration space GL + (3) and collect some results for 
the self-gravitating potential V. The potential V(F) at a typical configuration 
is an elliptic integral except for symmetric configurations (i.e. with at least two 
equal singular values) for which it can be integrated by elementary functions. 
The extension to the unconstrained space of the potential V depends on F e 
GL + (3) through the two principal invariants I\ and I2 defined by © and (HJ) 
respectively. In the study of the existence and stability of relative equilibria of 
simple mechanical systems will be necessary to compute the first and second 
derivatives of V^ L ^ a and the results of this section are essential to this end. In 
the following we will restrict ourselves to diagonal configurations only. There 
is no loss of generality with this assumption since, according to the singular 
value decomposition, every matrix F 6 GL + (3) belongs to the G-orbit of some 
diagonal configuration D by some element (e; A, B) £ G. Hence by Remark |3. II 
the qualitative properties of a relative equilibrium (D, £_r)) are the same as 
those of (ADB T , (A£ L ,B£ R )). 

Let J(k,r), with k,r € N be the following family of integrals: 



J(k,r) 



s r ds 
A(Fp 



and denote by V% (i = 1, 2) the partial derivative of V with respect to /j and 
by Vij the partial derivative of V% with respect to Ij for j — 1,2. Using ([6|), 
elementary calculus computations give: 

V = -RJ{1,Q) 7i = |j(3,2) y 2 = |j(3,l) (16) 

V n = -™J(5,4) V 13 = ~J(5 > 3) V 22 = -^J(5,2). (17) 

Note that the integrals J(k,r) are all positive as well as V\ and V%. 

Next proposition gives the values of J(k, r) in the case of spheroidal (two 
equal axes) configurations. 

Proposition 4.1. Let F — diag(a, a, c) be a spheroidal configuration with a 2 c = 
1. 

(i) The integral J(k,r) for the oblate spheroid F (a > c) and eccentricity 
e = \J 1 — ) 2 is given by 

9 f 1 (\ - r 2 Y T 3(fe-l)-2r J 

Mk,r) = 2 k / (1 (18) 
(1 - e 2 ) s Jo (l-e 2 x 2 )2 

(ii) The integral J(k,r) for the prolate spheroid F (a < c) and eccentricity 
e 



(f) 2 is given by 



rl (1 - x 2 ) r X^ k -V - 2r dx 

(l- e 2)^±^y (l-e 2 * 2 )* 



Jp{k,r) = - -—^ . (19) 
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Proof. Note that for the diagonal configuration F — diag(a, b, c) the value of 
A(F) in definition is 

A{F) = [( a 2 + s)(b 2 + s)(c 2 + s)}^ 2 . 

For (i): making the change of variables s — a? tan 2 9 we get 

Jo(k, r) = 2a 2 ^)-™ / (Smg) 2 (C 2 ° Sg) fc/2 d». 

Since the eccentricity is e 2 = a "^ c2 then a = (1 — e 2 ) -1 / 6 and c = (1 — e 2 ) 1 / 3 
because a 2 c = 1. Then, from the above expression for Jo(k,r) one gets 

r, /.Tr/2 / • /]\2r+l / m 3(fc-l)-2r 

(l_ e 2)2fc±i)^ 7 Q (l- e 2 cos 2 f?) fc/2 

Making x — cos 6* the result follows. 

For (ii): The change of variables s — c 2 tan 2 6 gives 

The eccentricity e of the prolate spheroid is such that c 2 (l — e 2 ) = a 2 . As 
a 2 c = 1 then a = (1 — e 2 ) 1 / 6 and c = (1 — e 2 ) -1 / 3 and the result follows for 
x = cos 0. □ 

As stated in the previous section, in order to find critical points of a G- 
invariant function defined in SL(3) we will work with its extension to GL + (3) 
subject to the constraint det F = 1. For, since any such function can be written 
as f(F) = /(7i(F), 12(F)), in order to compute the critical points we use the 
augmented function f x (F) = f(Ii(F), 12(F)) — X det (F) subject to the condition 
det(F) = 1. For the differentiation of / A consider the pairing between vectors 
SF G T F SL(3) and covectors B e T£SL(3) to be B ■ SF — tr(B T SF). Then 
using the chain rule we get that critical points must verify the following set of 
equations: 

df x (F) ■ SF = 2tr ^/ X F T + f 2 (hF T - F T FF T ) - ^ det(F)F"M SF (20) 

= 0. 
det(F) = 1 

since 

dh(F)-SF = 2tr (F T 5F) (21) 

dI 2 (F) ■ SF = 2 [h tr (F T 5F) - tr (F T FF T 5F)} (22) 

ddet(F) • SF = de^F^F^SF). (23) 

(see for instance [4] or [12]). 

Next proposition gives the expression for the locked inertia tensor. 



13 



Proposition 4.2. The locked inertia tensor for the G-action on T*GL + (3), at 
a configuration F G GL + (3), is defined by 



^ VlFF T + Cr F t F^~Tl F^F t ~U F T r\i J F 



•where T = j^p and ffi g so(3) for i = 1, 2 



(24) 



J-O ■ " ' ■ _ _ _ 

Under the isomorphism £/&e locked inertia tensor is also equivalent to: 



(^l^r)AF)(Vl,Vr))=T[^l Zr] 



i s -2det 
-2det (F)^" 1 i c 



— T 



»7l 
Vr 

'(25) 

where S — FF T , C = F T F and \a — tr(A)I — A (I denotes the identity matrix). 

Proof. By the locked inertia tensor definition in Proposition ^. 1[ the expression 
(jHJ for the infinitesimal generators of the G-action on GL + (3) and the definition 
(fT|) for the Riemannian metric, we have 

=« TlF - fTr,vIF - Frf^ » 



: Ttr 



TlF-F^) (jHF-Frm) 



Using the fact that and ?fi are skew-symmetric matrices and the cyclic prop- 
erty of the trace of a matrix, it is straightforward to obtain expression (|24[) . 

For the expression (|25p we need some standard properties of the isomorphism 
(J5J. In particular, 



tr 



tr (^77) =2e-ry 
-tr (1{L-I? 



14 = det(L)L- T ^X- 1 



(26) 
(27) 

if L is a symmetric matrix (28) 
if L is an invertible matrix (29) 



where • denotes the standard inner product on R 3 . Let us compute some terms 
of the expression (|24|) since the other are done similarly 



tr 



(zJrtiFF 1 



1 



tr 



Tl (vlFF T + FF T ij2) 



T - 



(by m) 

(by 
(by 
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tr (tfF^F) = det( ^ 1} tr [^F" V] (by (29}) 

= 2det(F)£ fl -F- 1 7 7i (by©) 



□ 

As a straightforward consequence we can obtain the momentum of a relative 
equilibrium for Dirichlct's model, that is its angular momentum and circulation. 

Corollary 4.1. The momentum for a relative equilibrium with configuration F 
and velocity-vorticity pair (£l>£r) £ R 3 ffi M 3 is 

M = I(F)(a,6?.) = T (i s a - 2det(F)F- T 6 ? , iota - 2det(F)F~ 1 a) ■ (30) 

XTiai is, the angular momentum and circulation of a Riemann ellipsoid with 
configuration by F, and angular velocity-vorticity pair (£l,£r) are given, re- 
spectively, by 

)/T = i s a-2det(F)F- T £ R 
c/T = i c ei?.-2det(F)F- 1 a 

The expression for the twice augmented potential ^ R follows now easily 
from Proposition 14.21 

-T {fa ■ isZL + fa ■ IcHr - 2 det(F)^ • F- T ^ R ) (31) 
-Adet(F). 

5 Existence conditions for symmetric Riemann 
ellipsoids 

In this section we classify symmetric relative equilibria for Dirichlet's model. We 
will treat the spherical case (i.e. a configuration having three equal principal 
axes) as a particular case of a symmetric configuration. From the singular value 
decomposition and the definition of the action of G = TL-i ix (SO (3) x SO (3)) 
on M = GL + (3) (or on SL + (3)) it follows that the stabilizer of a symmetric 
configuration F is conjugate to the stabilizer of a diagonal configuration. That 
is, conjugate to Z2 k 0(2)^ or Z2 x SO(3) D if F has 2 or 1 different singular 
values, respectively (see [20] for a derivation of this result). Actually if the 
configurations are diagonal, these are exactly their stabilizers. Here, K D denotes 
the diagonal embedding of K C SO (3) in SO (3) x SO (3) and 0(2) e is the 
subgroup of SO(3) generated by all the rotations Rg £ SO(2) e around a given 
axis e in M 3 and a rotation, II e _L , by tt around an axis e 1 - perpendicular to e. 
In case of the diagonal configuration F = diag(a, a, c), Rg is the rotation matrix 
by and angle 9 around (0, 0, 1) and Ii e ± can be chosen to be diag (1, — 1, — 1). 
We introduce the following subgroups: 
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• S0(2) e x S0(2) e , generated by elements (e; Rg 1 , Rg 2 ), with Re 12 € S0(2) e 
and (a;U e ±,Il e ±), 

• 0(2) e , generated by elements (e; Re, Re), with Rg g S0(2) e and (ex; U e ± , U e ± ), 

• Z 2 (e), the cyclic group isomorphic to Z 2 generated by the element (e; n e , n e ). 

• More generally, if if is a subgroup of SO(3) x SO(3), we denote also by K 
the subgroup of Z 2 x (SO(3) x SO(3)) generated by elements (e;k), with 

fee K. 

Note that since we are going to impose the constraint F € SL(3), we will 
consider only two kinds of symmetric configurations, especifically: 



spherical: F = diag(l, 1,1), 

spheroidal: F = diag(et, a, c), with a 2 c = 1. 



To find all the possible Riemann ellipsoids with symmetric configurations, 
we will have to solve (fT5|) with F of the above forms and different pairs (fa, fa). 
The possible solutions are summarized in the following theorem. 

Theorem 5.1. The relative equilibria for Dirichlet's model of a self- gravitating 
fluid mass are: 

(i) The spherical equilibrium with spherical configuration F — diag(l, 1, 1), 
velocity-vorticity pair (0,0) and Lagrange multiplier A = 2Vi + Its 
corresponding momentum and isotropy groups are 

H= (j,c) = (0,0), G M = Z 2 x(SO(3)xSO(3)), G F = G PF = 1 2 xSO(3) D . 

(ii) The family of MacLaurin spheroids which have oblate spheroidal confi- 
gurations F = diag(a, a, c) (with c < a) and angular velocity and vorticity 
parallel to the axis of symmetry . In terms of the parameter O defined 
by flea — fa — fa, this family is characterized by A = 2(1 — e 2 ) 2 l z V\ + 
4(1 — e 2 ) 1 ' 3 V2 and the following constraint between f2 and the eccentricity 

" 2 „VT^, ~2,„„. 6 n 2 . 



= 2- — (3-2e 2 )arcsine-— (1-e 2 ). (32) 



npG 

Rs corresponding momentum and isotropy groups are: 



(j,c) =2(l-e 2 )- 1 / 3 ™(e 3 ,-e 3 ), 



G F = Z 2 x 0(2)^ 3 , G,=SO(2) e3 xSO(2) e3 , G PF = 0(2) e 
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(Hi) Two branches of transversal spheroids which have prolate spheroidal 
configurations F — diag(a, a, c) (with c > a). We distinguish the two 
branches of this family with the signs + and These branches are char- 



acterized by the Lagrange multiplier A = 2((1 — e 2 ) 1 / 3 V A i + (l — e 2 ) 1 / 3 (e 
2)V2), the velocity-vorticity pair (£l,£r)± = w±(n, /±n) with f± 



2 _ 
l±e 



f where n is a unit vector perpendicular to e3 ) and the following constraints 
between w± and the eccentricity e: 

"1 _ (eTl) 2 (e±l) 



T ■> (3e + (e 2 - 3)arctanhe) . (33) 



irpG 

The corresponding momentum and isotropy groups are: 
( e(e±2) e (e T 2) 



(1-e 2 ) 2 / 3 ' (eTl)(l-e 2 ) 1 / 6 



G F = Z 2 k 0(2)£ , G p = SO(2) n x SO(2) n , G PF = Z 2 (n). 

Before proving the theorem we remark that formula (|32p has already been 
obtained by MacLaurin in 1742, as it is claimed in page 4 of Chandrasekhar's 
book [3]. 

Proof. First, using (J3J), (fTrJ|) and ([2U)l . it is easy to see that the general condi- 
tions ([T5)l are equivalent to 

= 2ViF T + 2V 2 {hF T ~-F T FF T )~\Aet{F)F~ 1 (34) 

-t [(Ul\\ 2 + Ur\\ 2 )f t - f t (£ l ® a) - ce« ® u)f t 

+2det(F) ({F-^l ® F~ t ^r) - (fr • F _t ^_r)F -1 )] 

1 = det(F) (35) 



spherical case: If F = I, then (|34l) . (|35| are simply 

(2Vi + 4V 2 - A — T [HClII 2 + H^ll 2 ~ ■ &] ) I+T Kl ® & + & ® & - 2^ <8> = 

Taking £l = £l,2, ££,3) an d the same sorte of notation for £r, the off- 

diagonal terms of this expression are independent of Vi , V2 and A, and equivalent 
to the following 6 equations: 





- £r,i)(£l,2 


- &,a) 


= 0, 


6i,l£z,,2 


= £rMl,i 




- 6j,i)(£l,3 




= 0, 


6i,l£z,,3 




(a, 2 


- £r,2)(£l,3 


- C-R.3) 


-0, 


£r,2£l,3 


= £r,3£l,2 



It follows then that £x, = £r- Recall from Theorem 12 . 1 1 that the momentum of a 
relative equilibrium with configuration x and velocity £ is given by /i = I(;r)(£). 
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Then, from we have fj, = (j,c) = (0,0). Therefore, G F = Z 2 k SO(3) d , 
G^ = G, G pjr = G^ n Gi? = G F . 

Now, noting that (£l,£r) £ flp F if £l = and that the velocity of a 
relative equilibrium is defined only up to addition of elements in g PF , the relative 
equilibrium (I, (£l,£r)) is the same as (I, (0, 0)). Then the relative equilibrium 
conditions are satisfied with A = 2 V± + W 2 • 

spheroidal case: We now consider F = diag(a, a, c) with a 2 c = 1. Since in 
the (ei,e2) plane all directions are equivalent, we can assume without loss of 
generality that = 0- Now, conditions (|3"4")l . (|3"5|) are equivalent to the system 



= 2aV l + 2a 3 V 2 + 2ac 2 V 2 - acX - aT(& 2 + &, 3 + £|, 2 + £|, 3 ) (36) 

+ 2T(c^, 2 a,2 + aUMs) 

= 2aVi + 2a 3 V 2 + 2ac 2 V 2 - acX - aT{? R1 + £| 3 + £f )3 ) (37) 

+ 2aT^, 3 a,3 

= 2cVi + 4a 2 cV 2 - a 2 A - cT(£ 2 ;1 + £|, 2 + £l, 2 ) + 2aTfo a &, 2 (38) 

= £r,i£r,2 = £ft,l£ft,3 = £ii,l£L,2 = £r,i£l,3 (39) 

= cT6^ 2 6w - 2aT6w£L,2 + aT& |2 &,3 (40) 

= aT£ fl , 2 £ fl)3 - 2aT^ )2 a, 3 + cT& )2 £ ii3 (41) 

1 = a 2 c (42) 



Note that these equations imply that = 0. Indeed, if ^ then (|39[) 
implies £r, 2 = £r, 3 = £l,2 = £ i:3 = and so (J2i])and ([37} imply = which 
is a contradiction. We will now proceed systematically considering four main 
cases: (i) £l = £r = 0, (ii) £l,2 = 0, (iii) £l,3 = and £ L ,2 7^ 0, and (iv) £ ii3 7^ 
and £l,2 7^ 0. 

(i) If £,l = £r = 0, then (f37|) . (|38l) and ([42]) are the only non-trivial conditions, 
which are equivalent to 

a 2 c=l, A = 4(2^1 + 4a 2 V2), (a 6 - l)(Vi + a 2 U 2 ) = 0. 

As V\ and are positive and < a 7^ 1 it follows from the last equation that 
there is no solution. 

(ii) If Ci,2 = 0, then from (36|) and (J37J) we have £ R>2 = and so ((38]) and (g2j) 
give A = 2c 2 V x + 4cV 2 . 

Let f2 = ££.3 — Then the remaining non-trivial equations, (|36p . (|3T[) 

and d42j) , give 

fi2= 2^-c^ 2 
i er 

As Vi and V 2 are positive, then last equality implies that the spheroidal config- 
uration i 7, is oblate, that is a > c. The eccentricity of the spheroid is e = 1 — ^ 
and a 2 — (1 — e 2 )" 1 / 3 . Using the relations $W\i for the partial derivatives of 
the self-gravitating potential, V% and V 2 , the expression Jq for the integrals 
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J(k,r) in the oblate case given by (fT5)l and & — 2npG, then the equation (|43[) 
is equivalent to 



|e 2 (j (3,2) + (l- e 2 )- 1 / 3 J (3,l) 
2 7 r / 9Ge 2 ( J (3, 2) + (1 - e 2 )" 1/3 J (3, 1 



One can easily compute the definite integrals Jo (3, 2) and Jo (3,1), although 
we avoid to display their expressions since they are quite lengthy. However the 
expression J (3, 2) + (1 - e 2 )' 1 / 3 J G (3, 1) is 



Jo(3, 2) + (1 - e 2 )- 1 / 3 J G (3, 1) = 3 ^-j— + — 7 — v 7 ! - e 2 arcsin e, 



2 



from which (|32|) follows. From Corollarv l4.1l it is trivial to obtain that the mo- 
mentum of this relative equilibrium. Hence, using (|9|), the appropriate isotropy 
groups are also straightforward. 

(iii) If £ Li 3 = and £ L , 2 + 0, it follows from flU and (gTJ) that £r >3 = 0. Then 
([3T| and 02]) give 

a 2 c=l, A = 2a 2 (Vi + (a 2 + c 2 )V 2 )- 

We can set n = e 2 and £r) = w(n, /n). So, substituting the above value of 
A into (l36|) and ([38]) these equations are 

= (a - 2c/ + af)Tuo 2 (44) 
= 2(c 2 - a 2 )^ + 2a 2 (c 2 - a 2 )V 2 - (c 2 / 2 - 2ac/ + c 2 )Tu? (45) 



From (|44|) we obtain the solutions /± = c±v '^ 2 a2 , from which follows that the 
spheroids must be prolate (c > a). In terms of the eccentricity e 2 = 1 — p- we 
have 

l±e 

J± = 



VT^ 2 

Therefore, the equation (|45|) gives 

2 2(c 2 -a 2 )(yi+a 2 V2) l^e 
W± ~ T( C 2/|-2ac/± + c 2 ) ~ T 



(Vt + (1 - e 2 ) 1 '^ 



Substituting in this expression V\ = yjp(3, 2) and V2 = -|jp(3,l), as well as 
R = 2npGT gives (|33|) . As before, using the expression of the locked inertia 
tensor and (|9|) the remaining results follow. 

(iv) In this case we have £l,2 7^ 0, £l,3 ^ and £z,,i = £^1 = 0. Note that from 
lHQj) and (|4*Tj) one should also have £^2 ^ and £^3 ^ 0. So, let £l,2 = g£,R,2 
and £2^3 = h^ji ^ for some reals g, h ^ 0. Then, using also (|42p . the equations 
([40)1 and gl]) reduce to 

(1 - 2a 3 5 + a '.,!, , = 

(a 3 - 2a 3 h + gh ;i /,..,;•/,■. ■, = 0. 
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These equations have solutions 



(h±,g±) 



'5a 6 - 1 ± Vl - 10a 6 + 9a 12 1 + 3a 6 ± \/l - 10a 6 + 9a 
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> 4a 6 ' 4a 3 I 

(46) 

Comparing (05]) and (J37J) using £ Li2 = g^R,2, £c,3 = ^6i,3 and (g2j! we obtain 
-a(l + g 2 ) + 2cg = and so g± must be 1+ ^l~ a6 or x ~^i~ g6 , but this is a 
contradiction with (|46|) . So there is no solution for the above system. □ 

Remark 5.1. Note that the functions /+ and /_ appearing in the charac- 
terization of transversal spheroids are inverse of each other. Therefore the 
two families of transversal spheroids belong to a single orbit of the symme- 
try group. Indeed, the Z 2 symmetry interchanges the + and — families, since 
u- {F\ w+n, uj_i_ /_|_n) = (F T ; w + / + n, us+n) = (F; w_n, w_/_n), as it follows from 
their definitions that a/j./w 2 = /_//+. 

Remark 5.2. 

1. Theorem 15. II is in agreement with Riemann's classification of ellipsoidal figu- 
res of equilibrium for Dirichlct's model of self-gravitating fluid masses. That is, 
these ellipsoidal figures of equilibrium must lie in one of the following categories: 
(a) the case of a uniform rotation with no internal motion (or uniform vorticity 
and no rotation); (b) the case when the directions of the angular velocity 
and vorticity £r are the same and coincide with a principal axis of the ellipsoid 
(also known as ellipsoids of type S); (c) The case when the angular velocity and 
vorticity are not parallel but lie in the same principal plane. 

In particular, we show that for Dirichlet's model it is not possible to obtain 
relative equilibria with spheroidal configurations belonging to category (c). 

2. The existence of transversal spheroids is referred in Chandrasekhar's book 
[3 a (see for instance page 143), however their study is not present in the classical 
works of Liapunov [TU] and Poincare [IS] • 



6 Stability conditions for symmetric Riemann 
ellipsoids 

In this section we apply the singular version of the reduced energy-momentum 
method introduced in [21] in order to deduce the stability of the symmetric 
relative equilibria obtained in Theorem 15.11 In order to apply this method it 
is essential to compute the second derivative of the twice augmented potential 
V c x c . Next lemma gives that result. 

Lemma 6.1. If F is a critical point of the twice augmented potential 
v Lin) = V{h{F),h{F)) - I [a Cr] 1(F) 



- Adet(F), 
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for (£l,60 € R 3 x R 3 and 1(F) as m Proposition gg| and A,Be T F GL + (3), 
then 



1 



df-V^^B) = d|y(A,B) - - [a (D|I(A,B)) 



where 



A det(F) (tr tr A) - tr {F^BF' 1 A)) 



d%V(A,B) = 2tr(B T A)(V 1 +I 1 V 2 ) 

- 2tr (BF T FA T + FB T FA T + FF T BA T ) V 2 

+ 4 tr (F T A) tr {FB T ) (V 2 + Vn + 2hV 12 + l\V 22 ) 

- 4tr (FF T FB T ) tr (F T A) (V 12 + J^a) 

- 4tr (F T FF T A) tr (FB T ) (V 12 + hV 22 ) 
+ 4 tr (F T FF T A) tr (FF T FB T )V 22 



and 



e«]D|.I(A,B) 



Ttr (4^5 T £^4 - AB T - 2^ B T A) 



Proof. We will just sketch the computation of d^V^-A, _B). 
Recall from the proof of Theorem 15.11 that 

dV(F) ■ 4 = 2(Vi + Z1V2) tr (F T A) - 2V 2 tr (F T FF T A). 

Differentiating again using the expressions for dli (F) ■ A and dl 2 (F) ■ A given 
in pTjl and ([22]) and the chain rule the result follows. 

For the expression T) 2 F 1(A, B) we differentiate the expression (|2~4"|) . which in 
this case takes the form 

I(*0(&,&)) =Ttr [2£> T £> - - FF^ 

Then, applying standard properties of the trace, we get 
(TlTr), (DI(F) • A)fa,bt)) - Ttr UTrA T TlF - 2A T fTr - 2AF t Tl 



The expression for D|J(A, B) stated follows now easily. Finally, using (f2"5| to 
differenciate the expression 

d det(F) ■ A = det(F) tr (F" 1 A) 



we get 



d% (det(F)) (A, B) = det(F) [tr (F~ X B) tr - tr {F^BF' 1 A)] . 



□ 
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6.1 Spherical equilibrium 



We now study the stability of the spherical equilibrium. Notice from the outline 
of the method in Section 3 that for this equilibrium whe have that q^, the cor- 
rection term and the Arnold form are all trivial, as well as the velocity- vorticity 
pair (£z,£r). As a consequence, sf,^ 3 " 1 = S SL ( 3 ), the orthogonal complement to 
the G-orbit at the identity in SL(3). Hence, to conclude stability of the spherical 
equilibrium we need to study the definiteness of 



dfV, 



(0,0) 



,SL(3)- 



Theorem 6.1. For Dirichlet's model, the spherical equilibrium is nonlinearly 
G-stable. 

Proof. Recall that TiSL(3) is the space of traceless matrices. Also, the infinites- 
imal action of q on GL + (3) at I is (£z,,6j)gl+(3)(I) = €l - £r- Then 

S SL < 3 > = {A G TiSL(3) :tr(4I-4^)=0 V£ L ,&ieR 3 }. 

Therefore S SL ( 3 ) is the space of traceless symmetric matrices. That is, matrices 
of the form 

an «12 a l3 
A = ai2 a 22 a 23 

_ai3 «23 -(aii + 022)_ 

We fix a basis for S SL ( 3 ) with respect to which the components of A are (an, a 22 , ai 2 , ai3, a 2 3) 



By Lemma RTTl the expression of dj Q ^ 
in Theorcml5.1i reduces to: 



gSL(3 



. , with A = 2 V\ + 4V 2 as given 



d\V, x 



(0,0) 



3) (A,B) =4(Vi + V 2 ) tr (BA), 



where we have applied the fact that A and B are traceless symmetric matrices. 
Therefore in this basis we have 



djV ( x 



(0,0) 



gSL(3) 



= 4(Vi + V 2 ) 



2 1 

1 2 













2 

2 

2 



As Vi + V2 is positive, the eigenvalues of this matrix are 12(Vi + V2), 4(Vi + V2) 
and 8(Vi + V 2 ) with multiplicities 1, 1 and 3 respectively. These are all positive, 
therefore the spherical equilibrium is G-stable. □ 



6.2 MacLaurin spheroids 

We now study the nonlinear stability of MacLaurin spheroids in the setup of 
previous sections (Theorem 15. ip . As it has been stated, a MacLaurin spheroid 
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has an oblate configuration which, with no loss of generality, we suppose diag- 
onal. This configuration is uniquely characterized by the eccentricity e G (0, 1). 
In order to apply Theorem 13. II one needs first to split g — R 3 ©R 3 according to 
(fTU|) . that is as 

Qn = 0PF © P and = &F © P © t. 

Recall that for the MacLaurin spheroid one has Gf = Z2 k 0(2)-° and G^ = 

SO(2)e 3 x SO(2)e 3 . One can then choose the following ordered orthonormal 
bases (with respect to the Euclidean product in R 3 © R 3 ) for each of the spaces 



of the above splitting: If we define h — -^{e s ,e 3 



P 



72 ^ 



-e 3 ),ti 



(ei,0), t 2 = (0,ei), t 3 = (e 2 ,0), U = (0,e 2 ), then 

q f = span{/i} 

p = spanjp} 

t = span{£i,i2,i3, U}- 

It is straightforward to check that these subspaces are invariant for Gp F — 

0(2)e 3 . The orthogonal velocity £ for the MacLaurin relative equilibrium is 
the orthogonal projection of the velocity £ onto p. Then, 

{£l,£r) X = ^(£L,3e 3 ,£fl,3e 3 ) ■ (e3,-e 3 )(e3,-e 3 ) 



?R,3)(e 3 , -e 3 J = 2"(e 3 , -e a ) = -j=p 

where f2 must satisfy (|32l) . As already defined, Iq is the restriction of I to p © t. 
The locked inertia matrix for the configuration F = diag(a, a, c) is, according 
to ((251), 

\D 1 D 2 



1(F) = T 



D 2 D 1 



where D\ — diag(a 2 + c 2 



- c 2 , 2a 2 ) and D 2 = — diag(2ac, 2ac, 2a 2 ). It is now 



straighforward to obtain the lo matrix with respect to the basis (p, t\,t 2 ,t3, ti). 
That is 



~2ac 
,2 1 ^,2 



= T 



4o 2 







a + c 
— 2ac 





a + c 





Or, in terms of the eccentricity 





2-e 2 



T 









.2^13^2 









-2ac 







— 2ac 
a 2 + c 2 



-2^r 







-2^r 

2-e : 



(47) 
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We can use Iq and (^l,£,r) ± to compute the momentum of a MacLaurin 
spheroid. Indeed 

- , n i ^ n , s 2V2Tfl 
A* = Wi>£«) = T/f W = (1 _ e2)1/3 P 

which is of course the same as the value obtained in Theorem 15.11 under the 
identification g — g* induced by the Euclidean product in M 3 © R 3 . 

In order to apply Theorem 13.11 we need to verify that the singular Arnold 
form is nondegenerate. 

Proposition 6.1. For a MacLaurin spheroid q M = t and the Arnold form, 
defined in (|14p , is positive definite for all eccentricities. 



Proof. Recall that the Arnold form Ar : q p x q M — > R is defined by: 
Ar( 7 i, 72 ) = (ad^/x, A(F,m)(72)>, 

where 

A(F,m)(7) =V (ad» + P p . ffii . [ad 7 (%VJ 

First, we compute the space q^. Notice the following relations for the adjoint 
representation of G: 

ad tl p = "y=^3, ad f2 p = ;4|*4, &d t3 p = ;4=*i> aA uP = ~/|*2> ( 48 ) 

and 

&d tl t 2 — ad tl i 4 = ad t2 t 3 = ad t3 i 4 = 0, 

ad tl t 3 =^(/i + p), (49) 

ad t2 i 4 = 7f( /l_ -P)- 

Also, recall that under our identification g ~ g* we have ad*p = — ad 7/ o, for 
7, G g, p G g*, and where p in the right hand side is identified with an element 
of g. Therefore P BF (ad* fi) = for i = 1, 2, 3, 4, hence q^ = t. As ^%{p) = M 

then Iq 1 (fj) = -^Lp. Then, from (|48j) we obtain 

ad tl (l V) = ~^*3 ad t2 (lo V) = i:U 

ad t3 (Iq V) = ad t4 (%V) = -7^2 

The inverse of the matrix (T47|) is not difficult to compute. Here we just state 
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the values of I 1 (ad^/i) = — 1 1 (ad^/Lt) on vectors w of the fixed basis: 

V 1 «M) = {1 _ e2 y /3 ^\h) = .4 *3 + ^4 *4 

-1/ „ \ 2TO -4vT^l?n 2(e 2 -2)f2 
lo 1 « M ) = - (1 _ e2)1/ 3 I o' 1 (^) = — ^ ta + g4 ^ 4 

- 1 , „ v 2TO - , , s 2(e 2 - 2)0 4\/l -e 2 ft 
V 1 «m) = - (^^ ^fti) = 1 e4 J *i - ta 

-1, \ 2TO , 4%/T - ^ 2(e 2 -2)ft 
\ ( ad ^) = (iT^jVa^ = * - e 4 t2 



From these expressions it follows easily that 



4yT 



e 



2 



A(F, M )(t!) = ^ — nt 3 + — nt 4 

A/rn w . -4\/T — i 2 ^ e 4 + 4e 2 -8^ 
A(f»(t 2 ) = — I-j nt 3 + 



4yT 



e 



2 



A(F )M )(t 3 ) = — ^ — mi - — nt 2 



. w s 4VT Tr ?^ e 4 + 4e 2 -8^ 

A(F,/x)(t 4 ) = — nti — 4 — nt 2 

Finally, the entries of the Arnold matrix are given by 

9 /OTO 

Ar(*i,*j) = (Mj adt 4 A(F , fx) (tj )) = _ /3 (p, ad f< A(F, M )fa)), t = 1,2,3,4. 
Using (I49[) the Arnold matrix is then given by 



Ar = 





-A 2 








-A 2 


At 














Ai 


-A 2 








-A 2 


Ai 



with 



A 1 = ^ nr 

e 4 (l-e 2 )TT 



A 8(l-e 2 )6TSl 2 
^2 — 7i 



The trace and the determinant of each block of Ar are positive and so Ar is 
positive definite. □ 

Next theorem gives the stability of the MacLaurin spheroids. 

Theorem 6.2. A MacLaurin spheroid with eccentricity e and momentum (i is: 

(i) nonlinearly G ^-stable if e < eo with eo — 0.952887. 

(ii) unstable if e > eo- 
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Proof. For (i): As the Arnold form is non-degenerate, G M -stability will follow 
whenever (d|,V^j. +corr^i(F)) is positive definite. For, recall from (|12p 
that 

Sint = {isl {3 )(F) + A : 7 £ q* 1 , i £ S SL < 3 \ and (DI(F) ■ ( 7sL(3) (F) + A)) (£ x ) G 

where S SL ^ 3 - ) is the linear slice for the G-action on SL(3) at the oblate configu- 
ration, F — diag(a, a, c). Matrices A G S SL ^ 3 ^ must verify 



= tv[A T ^F - Ff})], V&tjGK 8 , 
= tr (F~ 1 A) 



and 



because A must belong, respectively to the orthogonal complement to the tan- 
gent space to the group orbit through F and A G TpSL(3). These two conditions 
give that A must be of the form 



A = 



an O12 

Gtl2 a 2 2 

a 33 



with -(an + a.22) + -«33 = 0. 
a c 



Therefore we can describe S SL ( 3 ) as the set of matrices of the form 



.4 



ai + a 2 a 3 
a 3 ai — a 2 
-2^01 



(50) 



with 01, 02; 03 G K. Let the vector 7 G q M = t be 7 = (71, 72, 73, 74) with respect 
to the basis (ti,t 2 ,t 3 ,t4). Therefore, using we have 



7SL(3)(^)+^ 



oi + a 2 a 3 C73 - 074 

a 3 ai - a 2 — c 7 i + 072 

-073 + 074 071-072 ~2foi 



The set Ej nt is precisely the set of matrices Xsl(3)(F) + A for which 

(DI(F) • (7sL (3 )(^) + A)) (Cl^r) 1 - G p* = ( F + t)°, 

where (f) + t)° denotes the annihilator of f) + t. Using (Cl,^r) ± = S^P and 
it) G {/i, t\, t2, ta, £4} and differentiating (|24p , the computation of 

<«,, (DI(F) • (Ag L(3 )(F) (t L ,£+R) X ) 
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gives, in terms of the eccentricity e: 
(h, (DI(F) • ( 7SL(3) (^)+^)) (a,6i) X ) = 



(t u (DI(F) • ( JSU3) (F) + A)) (a, £fi) X ) = on ™ 2U/3 (2 + e 2 ) 73 - 2^1^74 



(t 2) (DI(F) • ( 7sL( 3)(F) + A)) fa,t a ) x ) = — ^ 2^1^73 - (2 + e 2 ) 74 



2(1 


_ e 2)l/3 




TO 


2(1 


_ e 2)l/3 




TO 


2(1 


_ e 2)l/3 




TO 



(t 3 , (DI(F) • ( 7S L(3)(^) + A)) (Zl,^) = 2(1 _ e2) i/3 [-( 2 + ^ + 2 V / T T ^72 
(t 4 , (DI(F) • ( 1SL{3) (F) + A)) ^ l ,Cr) X ) = 2(1 ™ 2)1/3 [-271^7! + (2 + e 2 ) 72 

It follows from the above expressions that 

(DI(F) • (A SL(3) (F) + A)) (a,^) X e (flF+t)° 

if and only if 7 = 0. This is equivalent to Ej nt = S SL ' 3 '. Let us now compute 
the correction term restricted to E int = S SL(3) . For any A G S SL ( 3 ) one has 
(DI(F) • A) (fc^fl) 1 - = ffffi oiP and so 

P t . fflp . [(DI(F) ■ A) (a, £*) X ] = ^ (1 4 ™^ /6 «i,0, 0, 0, 0^ . 

From the expression of lo it is straightforward to obtain 

% 1 (P t . er [(DI(F) • B) (fc, Zr) 1 -] ) = V2il(l - e 2 ) 1/6 biP, 
where b\ is the entry of B G S SL ( 3 ' playing the same role of a\ in A. Then, from 

CES 



corr (££ , ?fl) x(F)(A,B) = ^^^L aiP , ^£1(1 - e 2 ) x/ %^ = 8m 2 0l 6i. 

The computation of d 2 F V^ L ^^ ± (A,B) is lengthy but with no difficulties. 
Using Lemma \G. II we obtain 

(fe.^'^D^AB)) (Z L ,Z R ) ± )=ia 1 b 1 TSf 

d 2 F det (A, B) = -2c(3ai&i + o 2 & 2 + 03^3) 

We fix a basis for the slice S SL ^ 3 - ) in which the coordinates of A in ([50)1 are 
A = (01,02,03). With respect to this basis the matrix for di Vj„ is given by 

d|VL = diag(Di, D 2 ,D 2 ) with 

10 

—D x =o 4 (o 6 + 2)Vi + 3o 6 (o 6 - 1)V 2 + 4(o 6 - l) 2 V n + 8a 2 (a 6 - l) 2 V 12 
+ 4a A (a 6 -l) 2 V 22 , 

10 

— £ 2 =a ll Vi + a 6 (l-a 6 )V2. 
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Therefore, with respect to this basis we have 

(d|.t$L+corr € JL (F))| s = diag^x+GcA+era 2 , D 2 +2c\, D 2 +2c\) = diag(5i,5 2 , S 2 ). 

(51) 

For the MacLaurin spheroid we have, from Theorem l5.ll A = 2(c 2 Vi + 2cV 2 ) 
and Q 2 — ^e 2 (Vi + a 2 V 2 ). Then, in terms of the eccentricity, we have 



Si 



((3 - 4e 2 + e 4 )Vi + 3(1 - e 2 ) 2 / 3 V 2 + 2e 4 (l - e 2 ) 2 / 3 y n 
+4e 4 (l-e 2 ) 1 /-V 12 + 2e 4 T/ 22 ) 

((2 - 3e 2 + e 4 )l/i + (2 - 3e 2 )(l - e 2 ) 2 / 3 l/ 2 ) 



1-e 2 



Expressing the partial derivatives of V in terms of the integrals Jo(k,r) we 
obtain 

2R 



51 = — (9e(3 - be 2 + 2e 4 ) - \ / l-e 2 (27 - 36e 2 + 8e 4 ) arcsine 

5 2 = ((1 - e 2 )e(3 + 4e 2 ) - v / l-e 2 {3 + 2e 2 - 4e 4 ) arcsin 



The plots of Si and S 2 are shown in Figure [T] They show that Si is always 
positive in (0, 1) while S 2 has a root eo € (0, 1) being positive for e < eo and 
negative for e > eg. We used the Mathematica programing system for the 
numerical computation of eo to obtain eo ~ 0.952887. Therefore MacLaurin 
spheroids are G^-stable for e < e . 




Figure 1: The plots of Si (left) and S 2 (right) in R units. 



For (ii): When e e [eo,l) we study the spectrum of the linearized Haniil- 



tonian vector field Ly 



^ 2 



(d 2 /i£-L L.), where in our case N and u>n are as in 



Theorem 13.21 Also, from Proposition 6.2 in [21] we have 



(dlV)U = 



Ar 

( d F^(i >5R) ^ + corr (a,e fl )^)|s SL ( 3 ) 

< • > S SL(3)» 
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Let us start by computing the block -C •, • ^ s sl<3). of the above matrix, where <C 
■, • ^s SL ( 3 )* is t ne inner product on the dual of S SL ( 3 ) induced by the Riemannian 
metric on S SL ^ 3 ^. For, let as before fix the ordered basis {si, S2, S3} on the slice 
S SL ( 3 ) where 





"1 










"1 





0" 




"0 


1 


0" 


Sl = 





f 





S2 = 





-f 





S3 = 


1 
















-2Vl -e 2 

























Recalling that <§; A, B ^ s sl(3>= Ttr(A T B) then the matrix that represents 
<C •, • 3>s SL < 3 ) m the fixed basis is 



111 = 2T 



'3 - 2e 2 







1 
1 



Let {s^jS^jS^} be the dual basis of {si,S2,S3} under the identification of 
S SL < 3 )* with S SL < 3 ) using the inner product <C •,• 3>s SL < 3 )- I n this basis the 
induced inner product <C •, • ^s SL < 3 )* i s represented by TZ^ 1 . 



Let K 2 = ^V^ L (R)± +coTT { ^ R y(F))\^ t , then 
( d2 F h ({L,uy)\ 



N 



Ar 

n 2 



K1 1 



To compute Lh in the basis {t\, t 2 , £3, £4, si, S2, S3, s*, S2, S3} for N we use the 
formula for a; at given in Theorem 13.21 Let us now compute each of the blocks 

Of U>N- 

Recall that q p = t for a MacLaurin spheroid. Then, from Theorem 13. 2\ for 
7i , 72 G t we have: 



£(71,72) = ' ad 7l7 2 = 
Using (I49p . the matrix 3 is given by 

2TVL 



(1-e 2 ) 1 / 3 



p,ad 7l 72 



(l_e 2 )V3 



-f 



-f 
f 





Since for a MacLaurin spheroid Si n t = S SL ( 3 ) it is immediate from the definition 
of ^ in Theorem 13.21 that ^ is the zero matrix. 

We now compute the Coriolis term —dx^ L '^ R ^ L ■ For that we will obtain 
a concrete expression for the right hand side of the equality 

dx iU,U^ (X; Y ) = X{ X ^* )X (Y)) - Y( X ^'^ (X)) - X ^' M± ([X, Y]), 

(52) 
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with X,Y e X(GL+(3)). 

Start by considering for U, V £ 7/GL + (3) = L(3) the corresponding left- 
invariant extensions Xu,Xy € X(GL + (3)). We have Xu{F) = FU for every 
F 6 GL+(3). Recall that (Cl,^)^ = f(e 3 ,-e 3 ). Then according to the 
definition given in Theorem 13.21 we have 

x («,,«h)^ (Xu)(F) = <f tr ((ejF + Fei) T FU) 

(53) 

= tr (e^(F T FU + FUF T )) . 
It is straightforward to obtain 



X v ( X ^^ X (Xu))(F) = 

= _£M1 tr (eJ((F^) T F?7 + F T FVU + FVU F T + FU{FV) T )) (54) 
• tr (eS(V T F T FU + F T FVU + FVU F T + FUV T F T )) 



nr 

2 



Also, since Xjj,Xv are left-invariant vector fields, the identity [Xu,Xy] 
Xuv-vu holds and we have, from (I5S 



X 



([Xu,X v ])(F) = 



In order to compute — dx^ L '^ R ' 



tr (e^(F T F{UV - VU) + F{UV - VU)F T )) . 

(55) 



„ let F = diag(a, a, c) and A, B € Si n t- 
The unique left-invariant vector fields extending A and _B are Xp-i^ and Xp-i b 
respectively. Then, using (fB"2"]) together with ([54")) and ([53)1 it is immediate to 



obtain — dv^ 1 "^) I =0. Therefore, from Theorcm l3.2l thc symplectic matrix 

int 

wjv and its inverse are 










0" 




"s- 1 


" 


LU N = 








1 


<V = 





-1 







-1 







1 






So the linearized vector field is 



'H _1 Ar 

ft 2 








where in our basis 7^2 = diag(,5i, S2, £2) is given in (|5T|) . 

The block S _1 Ar has imaginary eigenvalues = ±i ^/SEiE Q with multiplic- 
ity 2. For the block 



-72; 



72a 







we obtain the following eigenvalues: 

(6-4e 2 )T w ith multiplicity 1, and 



5$ = ±i 
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= ±i\/j^ with multiplicity 2. 



As (6 — 4e ) and 5i are positive for < e < 1, then ef is always imaginary 
in that range. However, becomes real if S2 becomes negative, which happens 
exactly at eo as seen from the previous stability analysis. Hence, if e > eo the 
MacLaurin spheroid becomes linearly unstable, therefore unstable. This loss of 
stability corresponds exactly to a collision at of for e = eo, which passes 
from being pure imaginary to be real. □ 

Remark 6.1. 

1. In Riemann's work [19] some conclusions were made concerning the stability 
of Maclaurin spheroids by studying the existence of a minimum of a certain 
function G. The existence of this minimum was not done by studying the 
second variation of G. Riemann even says in page 188: "The direct analysis 
of the second variation of G when the first variation vanishes would be very 
complicated; we can however decide if this function has a minimum by the 
following form:...". He follows with the analysis of the behavior of G. His 
final conclusion on the stability ends the paragraph 9 of his paper and is the 
following: "From this study it follows that the case of a rotation of an oblate 
ellipsoid, around its shortest axis, case already known to MacLaurin, can only 
be unstable if the relation between the shortest axis with the others is less 
than 0.303327...". We note that if the relation between the shortest axis of the 
Maclaurin spheroid and the others is | < 0.303327 this is equivalent to say that 
the eccentricity e > eo = 0.952887. 

The value 0.303327 obtained by Riemann follows from his study on the 
existence of oblate spheroids in pages 184-185 of [19], namely as the root of the 
last displayed equation in page 184. We remark that this equation is equivalent 
to the equation S2 — where S2 is as in the proof of (i) in Theorem 16. 21 Indeed 



S 2 = <S=> e(l - e 2 )(3 + 4e 2 ) - yj\ - e 2 (3 + 2e 2 - 4e 4 ) arcsine = 0. 
Taking e = cos-0 = sin(~ — ip) and ip e (0, 5) one has 

/5 1 \ 7T 

cos ip sin 2 ip (5 + 2 cos(2-0)) — sin ip ( - — cos(20>) — - cos(4?/>) J (— — ^) = 
10 sin(27/>) + 2 sin(4-0) + (-5 + 2 cos(20) + cos(40)) (tt - 2-0) = 0. 

(56) 

The equation (|56"|) is the same one appearing in Riemann's paper. With respect 
to the computation of the root of this equation Riemann just says: "...this 
equation has, for ip between and 7r/2, the unique root sinip = 0.303327...". 
Indeed, this is equivalent to say that eo = cos ip = 0.952887... 
2. Concerning the stability of MacLaurin spheroids, it is shown in Liapunov's 
work [10] that under the hypothesis of the preservation of the ellipsoidal shape 
(setup we used) the value eo of loss of stability is the same as in Riemann's 
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work, but if this hypothesis is dropped then the MacLaurin spheroid is only 
stable for e < e\ with e\ — 0.8126.... The point e\ is exactly the point where 
the family of MacLaurin spheroids bifurcates into a branch of ellipsoids with 3 
distinct axes lengths (Jacobi ellipsoids). We refer the reader to pages 52 and 
61-63 of [TD]. 



6.3 Transversal spheroids 



All the qualitative properties, including stability, of two relative equilibria lying 
in the same orbit of the symmetry group are the same. Therefore, in view of 
Remark 1 5. 1[ in this subsection we analyze the nonlinear stability of the + family 
of transversal spheroids and the main result, Theorem 16.31 will follow for both 
families. We will set uj+ = lo, f = /+ and = (£l,£r) for notational 

simplicity. Also, to keep the notation consistent with the proof of Theorem l5.1( 
we will set n = e 2 . 

In this case the computation of the splitting (fT0|) is simplified due to the 
fact that g Px = {0} and therefore p = g M . Introducing the vectors h — 
e 2 ,0), p 2 = (0,e 2 ), ti = (exjO), t 2 = (0,e x ), t 3 = -L(e 3 ,-e 3 ), 



75(e 3 ,e 3 J, Pi 
we choose 



qf = span{/i} 
p = span{pi,p 2 } 
t = span{*i,t 2 ,*3}- 

These subspaces are obviously invariant under the action of Gp F = Z 2 (e2). 
With respect to the basis (pi,p 2 , ti,t 2 ,t 3 ) for p © t we have, 



In =T 



a + c 
—2ac 






—2ac 
a 2 + c 2 









a 2 + c 2 
—2ac 







— 2ac 
a 2 + c 2 











4a 2 



Or in terms of the eccentricity of a prolate spheroid, 







2-e 2 


-2 




(i 


_ e 2)2/3 


(l-e 2 )V6 






-2 


2-e 2 




(i 


_ e 2)Ves 


(!_ e 2)2/3 


T 





















































2-e 2 




-2 


(1 


_ e 2)2/3 

-2 


(1 


_ e 2)l/6 

2-e 2 


(1 


_ e 2)l/6 




(1 


_ e 2 J2/3 





4 (l_ e 2)l/3 _ 



It follows immediately that p and t are orthogonal with respect to lo, so our 
choice of the splitting g — 0F©P©tis correct. In this basis the orthogonal 
velocity is 

with lo and / as in the + family in Theorem 15. II 

It is straightforward to compute the adjoint representation of g in this basis: 
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Lemma 6.2. The elements of p © t satisfy the following relations: 
ad tl i 2 = ad tl i 3 = -^Pi &d t2 t 3 = -^p 2 
ad tl pi = —{h + 1 3 ) ad tl p 2 = 

adt 2 pi = ad t2 p 2 = -^(h - t 3 ) 

ad t3 pi = ~^7^ tl adi 3 p 2 = 7^2 



ad Pl p 2 = 0. 

With this, we can prove the following proposition, which shows that the 
stability method is applicable. 

Proposition 6.2. The Arnold form for a transversal spheroid is positive defi- 
nite. 

Proof. Recall that the momentum of a transversal spheroid is 

M = %(£l, £r)~ L = Tcj (k l Pi + K R p 2 ) , 

where k l := [(a 2 + c 2 ) - (2/a)/] and k r := [(a 2 + c 2 )f - (2/o)] . From Lemma 
16.21 and recalling that ad*p = — ad 7 p, we have, for 7 = 7^ £1+7^ £2+7^^3 G t: 

P flF (ad» = ^( 7 «« i +7 (2) « fl ) h. 
This is zero iff 7^) = — K^ l \ with 

(e-l)(e + 2) 



k := 



fl (e-2)x/T _ ? 



and therefore = {7 (1) (*i - +7 (3) ^3 : (7 (1) ,7 (3) ) e R 2 }. A basis for 
is given by {^ a — t\ — Kt 2 ,^fb = £3}- Now, proceeding as for the MacLaurin 
spheroid, we compute 

ad 7a (l V) = ^= ((1 - Kf)h + (1 + nf)t 3 ) ad 7b (lo V) - ^(/«a - *i) 
ad 7o /i = -\Z2rwK L t 3 ad^/j = 7= (^ii - k%) 

(ad! u) = — ~ K -i 3 

V 272(1 -e 2 ) 1 / 3 



% V» = ((-VT^e 2 - 2)^ + 2(e 2 - l)««) fl 

©(a/T - ^(e 2 - 2)k r - 2(e 2 - l)/t L )t 2 ) . 
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From where it easily follows that, in the basis { 7a ,7b} for q^, 

/ 3e 4 (2 + e)Tcj 2 4(1 + e)(2 - e)(2 + e)Toj 2 
r_ iag l v 2(2-e)(l-e 2 ) 5 /3' e 2 (1 _ e 2)2/3 

The entries of Ar are obviously positive. □ 

We can therefore apply the singular reduced energy-momentum method to 
study the stability of the transversal spheroid, obtaining the following 

Theorem 6.3. Both families of transversal spheroids are nonlinearly stable for 
all eccentricities. 

Proof. We start by computing the space of internal variations S; n t . Recall that 
the slice S SL(3) at F = diag(a, a,c) is given by l[50]). Hence, for 7 = (7 (1) ,7 (3) ) 
in q'' and A — (ai, a 2 , 03) in S SL ^ 3 - ) we have 



7sl(3)(-F) 



at + a2 »3 — 
03 + v / 2a7 < - 3 - ) a\ — a 2 — (c + an)^ 1 ^ 







(a + ck)^ 1 ^ 



-2c 



«1 



Differentiating (|2~4")) we find 



(h, (DI(F) • (7sL( 3 )(i ;i ) + A)) fe,^) =0 

(ti, (DI(F) • ( 7SL(3) (F) + A)) (fr.fr)-^ = Ve+iw^+^i-WV") 

(t 2 , (m(F) . ( 7sl(3) (f) + a)) ( & ,fc)-L> = - g "^3;-f' 1/V3 '' 

(t 8) (DI(F) • ( 7SL(3 )(F) + A)) (fc,^) = (f^f^a - 

In order for (DI(F) • ( 7s sl(3) + A)) {(l,^r) ± G (flF © t)° all the above expres- 
sions must vanish, which happens if and only if 7W = and 7 ' 3 -' = 036, with 
e = 7=; ^ 6 om /« ■ Therefore we can chose a basis for the space of internal varia- 

tions Sjnt such that any element v belonging to it has components (at, a 2 , 0,3) 
with the parametrization 



v = (01,02,03) 



(1 - Viae) 
a 3 (l + V2ae) ai - a 2 







-2<; 



«1 



It is straightforward to obtain 
(m(F)-v) (£ L ,Z R ) ± =2eTu J 



((e-l)ai + (e+l)o 2 ) m , ( e + l)ai + a 2 (e - 1 ) 

(e-l)(l-e 2 )V 8 



(1 - e 2 ) 5 /6 



-Pi 



-% 



Then the correction term, for V\ = {a\, a 2 , a 3 ), v 2 = (pi, b 2 , 63), is given by: 



(vi,v 2 ) 



8Tlo 2 {9-5e< 



1 



1 



oi&i - 3(ai6 2 + a 2 6i) - a 2 b 2 
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For the restriction of the Hessian of V% L ^ ± at F we compute first the following 
second variations 

d p V(vi,V2) = Didibi + D 2 a 2 b 2 + D 3 a 3 b 3 , with 

4 



D 3 



(-(1 - e 2 ) 2 / 3 (-3 + e 2 )V x + 3e 2 (e 2 - l)V 2 + Ae 4 V n 



(1 - e 2 ) 5 /3 
+8e 4 (l - e 2 )^V 12 + 4e 4 (l - e 2 f'^V 22 ) 



4 [V 1 



(1 - e 2 ) 2 / 3 



T ((e 4 -l)V4+e 2 (l-e 2 ) 1 / 3 (e 2 -3)V2) 



((a,60MD 2 F i(^2))(a,£fln 



4T^ 2 



( e -?L)^(i+ e ) a i^i + T^( a i fe 2 + a 2 bi) + jjz^a 2 b 2 + (1 + e)a 3 b 3 ^j 



dpdet (vx,v 2 ) = 



-2 



yyj (3ai&i + a 2 b 2 + (1 - e 2 )a 3 & 3 ) 



(l_ e 2)l/3 

Putting all the contributions together, and substituting the integral expres- 
sions for the derivatives of the potential and the values of A and u> given in The- 
orem l5.ll we find that the matrix representing d? F Vh £ R s± +corr^ i ^ R ^± (F)\^ 
in the fixed basis of E; nt is block diagonal, with a 2 x 2 symmetric block in the 
first two components, and the other block given by the coefficient of 0363, say 
<f>, which has the following expression: 



(1 - e 2 ) 2 / 3 (3 + e 2 )V x + (3 + 2e 2 - e 4 )V 2 



(1 - e 2 ) 2 / 3 



It is clear that <j> is positive since (3 + 2e 2 — e 4 ) is positive in (0,1) and V\ and 
V 2 are always positive. 



The 2x2 block is U 



x y 
V z 



with 



x = -27e(l - e 2 ) + (27 - 36e 2 + 17e 4 )arctanhe 

y = -9(1 - e 2 ) [3e - (3 - e 2 )arctanhe] 

z = -2e(3 - 4e 2 ) + 2(3 - 5e 2 + 2e 4 )arctanhe 

The matrix U has real eigenvalues. The trace and the determinant of U are: 
R 



tv(U) 



[-e(33 - 35e 2 ) + (33 - 46e 2 + 21e 4 )arctanhe] 



35 



and 

det(U) = ^ [-27e 2 (21 - 40e 2 + 19e 4 ) 

-2(-567e + 1269e 3 - 831e 5 + 121e 7 )arctanhe 

- (567 - 1458e 2 + 1212e 4 - 334e 6 + 13e 8 )arctanh 2 e] . 

These two quantities are positive in the interval e G (0,1). Their plots are 
displayed in Figure O Therefore, both families of transversal spheroids are 
stable for all eccentricities. □ 




Figure 2: The plots of tr(J7) (left) and e 10 det(C/) (right), in R and R 2 units 
respectively. 



References 

[1] Arms, J.M., Marsden, J.E., and Moncrief, V. [1975], Symmetry and bifur- 
cations of momentum mappings. Comm. Math. Phys, 78 455-478. 

[2] Arnold, V.I. [1966], Sur la geometrie differentielle des groupes de Lie de 
dimension inhnie et ses applications a Phydrodynamique des fluids parfaits, 
Ann. Inst. Fourier, Grenoble, 16, 319-361. 

[3] Chandrasekhar, S. [1987], Ellipsoidal figures of equilibrium, Dover Pub. 
Inc., New York. 

[4] Ciarlet P.G. [1993], Mathematical Elasticity, volume 1: Three dimensional 
elasticity, North-Holland. 

[5] Cohen, H., and Muncaster, R.G. [1988], The theory of pseudo-rigid bodies, 
New York: Springer- Verlag. 

[6] Fasso, F. and Lewis, D. [2001], Stability properties of Riemann ellipsoids, 
Arch. Rational Mech. Anal. 158, 259-292. 



36 



[7] Lerman, E. and Singer, S.F. [1998], Stability and persistence of relative 
equilibria at singular values of the moment map, Nonlinearity, 11, 1637- 
1649. 

[8] Lewis, D. Marsden, J., and Ratiu, T.S. [1987], Stability and bifurcation of 
a rotating planar liquid drop, J. Math. Phys., 28, 10, 2508-2515. 

[9] Lewis, D., Simo, J. C. [1990], Nonlinear stability of rotating pseudo-rigid 
bodies, Proc. Roy. Soc. London Ser. A, Vol. 427, 1873, 281-319. 

[10] Liapounoff, M.A. [1904], Sur la stabilite des figures ellipsoidales d'equilibrie 
d'un liquide anime d'un movement de rotation, Annales dc la Faculte des 
Sciences de L'Universite de Toulouse. 

[11] Marsden, J.E. [1992], Lectures on Mechanics, Lecture Note Series 174, 
LMS, Cambridge University Press. 

[12] Marsden, J.E. and Hughes, T.J.R. [1983], Mathematical Foundations of 
Elasticity. Prentice-Hall. 

[13] Montaldi J. [1997], Persistence and stability of relative equilibria, Nonlin- 
earity, 10, 449-466. 

[14] Ortega, J-P. and Ratiu, T.S. [1999], Stability of Hamiltonian relative equi- 
libria, Nonlinearity, 12, 693-720. 

[15] Patrick, G.W. [1992], Relative equilibria in Hamiltonian systems: the dy- 
namics interpretation of nonlinear stability on the reduced phase space, J. 
Geom. Phys., 9, 111-119. 

[16] Patrick, G.W., Roberts, M. and Whiff, C. [2004], Stability of Poisson and 
Hamiltonian relative equilibria by energy methods, Arch. Rational Mech. 
Anal, 174, 301-344. 

[17] Pcrlmuttcr, M., Rodrfguez-Olmos, M. and Sousa-Dias, M.E. The symplec- 
tic normal space of a cotangent-lifted action, Differential Geometry and 
Applications, to appear. Preprint math.SG/051207. 

[18] Poincare, H. [1885], Sur L'Equilibre d'une Masse Fluide Animee d'un Mou- 
vement de Rotation, Acta Mathematica, 7, 259-380. 

[19] Riemann, B. [1861], Ein Beitrag zu den Untersuchungen iiber die Bewegung 
cincs fliissigen gleichartigen Ellipsoides, Abh. d. Konigl. Gesell. der Wis. 
zu Gottingen, 9, 168-197. 

[20] Roberts, M. and Sousa-Dias, M.E. [1999], Symmetries of Riemann ellip- 
soids, Resenhas- IME. USP, Vol. 4, 2, 183-221. 

[21] Rodrfguez-Olmos, M. [2006], Stability of relative equilibria in simple me- 
chanical systems, Nonlinearity, 19, no. 4 , 853-877. 



37 



[22] Rodrfguez-Olmos, M. and Sousa-Dias, M.E. [2002], Symmetries of relative 
equilibria for simple mechanical systems. SPT 2002: Symmetry and per- 
turbation theory (Cala Gonone), 221-230, World Sci. Publ. 

[23] Rosensteel, G. [1988], Rapidly rotating nuclei as Riemann ellipsoids, Ann. 
Physics, 186, 230-291. 

[24] Rosensteel, G. [1988], Geometric quantization of Riemann ellipsoids, Group 
theoretical methods in physics (Varna, 1987). Lecture Notes in Phys, 313, 
253-260. 

[25] Rosensteel, G. [2001], Gauge theory of Riemann ellipsoids, J. Phys. A, 34, 
13, L169-L178. 

[26] Simo J.C., Lewis, D. and Marsden J.E. [1991], Stability of relative equilib- 
ria. Part I: The reduced energy-momentum method, Arch. Rational Mech. 
Anal, 115, 15-59. 



38 



